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ABSTRACT 

We report the discovery of a circumbinary disk around the Herbig Ae/Be system v892 Tau. Our 
detailed mid-infrared images were made using segment-tilting interferometry on the Keck-1 Telescope 
and reveal an asymmetric disk inclined at ~60° with an inner hole diameter of 250 mas (35 AU), 
approximately 5x larger than the apparent separation of the binary components. In addition, we 
report a new measurement along the binary orbit using near-infrared Keck aperture masking, allowing 
a crude estimate of orbital parameters and the system mass for the first time. The size of the inner 
hole appears to be consistent with the minimum size prediction from tidal truncation theory, bearing 
a resemblance to the recently unmasked binary CoKu Tau/4. Our results have motivated a re-analysis 
of the system spectral energy distribution, concluding the luminosity of this system has been severely 
underestimated. With further study and monitoring, v892 Tau should prove a powerful testing ground 
for both predictions of dynamical models for disk-star interactions in young systems with gas-rich disks 
and for calibrations of pre-main-sequence tracks for intermediate- mass stars. 

Subject headings: techniques: inter ferometric, stars: pre-main sequence, stars: binaries: close, in- 
frared:stars, stars: individual (v892 Tau) 



1. INTRODUCTION 

v892 Tau (Elias-1, Elias 3-1) is a young stellar object 
in the Taurus - Auri^ae star forming region {d = 140 pc, 
iKenyon et al.l il994) . Its visible spectrum is classified 
as spectral type B8V (jHernandez et al.]l2004f ). making 
v892 Tau one of the few Herbig Ae/Be stars in Taurus 
(in addition to AB Aur, MWC 480). The spectral energy 
distribution (SED) of v892 Tau is dominated by bright 
thermal emission in the mid-infrared ( Hillenbrand et al. 
[1992) , with relatively large line-of-sigh t extinction in the 
visible (one estimate is Ay ~ 5.9, iKenvon fc HartmannI 
[1995). Despite being one of the closest Herbig stars and 
after decades of mult i- wavelength observations, the na- 
ture of v892 Tau is still uncertain. 

The SED suggests v 892 Tau is either a young em- 
bedded Class I source (|Lada] [19871 ) still surrounded by 
its nascent envelope or a more evolved. Class II object 
(i.e, Herbig Ae/Be star) seen through its edge-on disk. 
Spatially-resolved imaging can easily distinguish between 
these two scenarios, and early speckle interferometry 
([Kataza fc Maiharal ll99TI : iHaas et al.lll997l : iLeinert et al] 
[2001) suggested the presence of an extended and elon- 
gated nebula in the near-infrared that could be due to 
scattering in bipolar lobes. 

The high-resolution speckle imaging of ISmith et al.l 
(|2005l ) clearly resolved the K band (A = 2.2/im) emis- 
sion to be coming from two unresolved stars in v892 Tau 
with little or no sign of a nebula or extended emission. 
The two stars of roughly equal brightness had an appar- 
ent separation of 55 milliarcseconds (7.7 AU) and some 
evidence of orbital motion was seen between epochs sep- 
arated by 7 years. 
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Although this binary should have a dramatic effect 
on the infrared emission, carving out a large hole in 
the circumbinary disk, recent ana lysis of the SED of 
v892 Tau including ISO data (Ack e fc van den Anckerl 
[2004) uncovered no distinct signature of the underlying 
binary^. Other workers (Liu et al. 2005. 2007) made use 
of the technique of nulling interferometry in the mid-IR 
(A = 10.3/im) to marginally resolve v892 Tau (FWHM 
~20 AU) along PA 164° consistent with normal (single- 
star) disk emission. 

Here, we report new mid-ilR imaging of the v892 Tau 
system which resolves these mysteries, discovering very 
extended and resolved emission that we interpret to be 
coming from a circumbinary disk. We also confirm the 
presence of the binary at K band and our new data allow 
first crude estimates of the orbital elements. Lastly, we 
report a new SED analysis which provides an improved 
determination of the system luminosity and our viewing 
geometry. 

2. OBSERVATIONS 

We observed v892 Tau using the Keck-1 telescope as 
part of two separate experiments, the Keck segment- 
tilting experiment and the Keck aperture masking ex- 
periment. Here, we briefly describe these experiments 
and give pertinent observing details. 

v892 Tau was observed on UT2004Aug31, UT 
2004Sep01, and UT2005Febl9 at the Keck-1 tele- 
scope using the 10.7/im filter (Aq = 10.7/im, 
A A = 1.55/im) of the Long; \y aveleng[:th Spectrometer 
(LWS, ICampbeh fc JonesI [2004D just before this instru- 
ment was decommissioned in 2005. In order to op- 
timize our calibration against changes in seeing (e.g., 
iTuthill e^^ l. 2000), we used the Keek's segmented pri- 
mary mirror in a novel "segment-tilting" mode, whereby 
we controlled the tilts and pistons of the individual mir- 



6 lAcke fc van den Ancke^ (|2004l ) did note strong 11.2/im PAH 
emission and anomalous infrared colors. 
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ror segments to create a set of 4 independent and non- 
redundant interference patterns (using 6 segments each) 
on the LWS detector focal plane. Furthermore, short ex- 
posures were used (tint = 90 ms< to, where to is the at- 
mospheric coherence time ~250 ms in the mid-iR) to ef- 
fectively freeze the atmospheric turbulence. For calibra- 
tion we interleaved target observations with calibrators 
a TdivJ , L Aur, and a Get. Details on this experiment 
and the implementat ion at Keck were first introduced by 
iMonnier e t al.' (2004) and more information ca n be found 
in recently published science papers (Weiner et al.l l2QQ6l : 
[Ireland et al.ll2QQ7l : [RaTagopal et al.ll2QQ7f ). 

The combined Fourier coverage of the four patterns 
over the three nights (total of 6 independent pointings) 
can be found in Figure [T] along with the visibility re- 
sults. The visibility-squared data and the closure phases 
were compiled and saved using the OIFITS data format 
(|Pauls et al.ll2QQ5i) and are available upon request. Fig- 
ure [1] shows the clear sign of highly-elongated disk emis- 
sion. 

We also obtained diffraction-limited observations of 
v892 Tau in the near-IR on UT2004Sep04 using the K- 
ban d filter (A =2.21/i m, AA =0.43/im) of the NIRC cam- 
era ([Matthews et al.l ll996) in conjunction with an annu- 
lus aperture mask placed in front of the secondary mir- 
ror. This observing mode was extensively utilized be- 
tween 1997 and 2005 (see most recent science papers, 
IMonnier et all l2007l : iTuthill et all I2008D and details of 
th e experimenta l desig n and performance can be found 
in 'Tuthill et ajj (|2000l ). We interleaved observations of 
the target with the unresolved calibrator 54 Per. The 
Fourier coverage and visibility-squared results are shown 
in Figure [TJ where the the binary nature is clearly re- 
vealed. 

3. ANALYSIS 

3.1. Image Reconstructions 

We used the BSMEM i mage reconstruction software 
(jBuscherl 119941 : iLawson e t al. 2004, 2006) for aperture 
synthesis imaging. Figure [2] shows the reconstructed im- 
ages for the mid-IR and near-IR data. In order to con- 
firm the asymmetric features in the mid-IR image were 
not due to artifacts in the BSMEM algorithm, we also 
used the independe nt image reconstruction code MACIM 
([Ireland et al.ll200 6) and found good consistency between 
the basic morphology, size scale, and level of emission 
asymmetry. As an additional data quality check, sep- 
arate images were made for each independent night of 
data and the resulting images all closely resembled the 
result shown in Figure O Note that we do not know the 
relative position of the near-IR image compared to mid- 
IR image and have presented each image centered in its 
corresponding frame. 

3.1.1. Mid-infrared: Elongated disk structure 

The mid-IR image shows elongated emission, approxi- 
mately 320 mas by 180 mas in full extent, elongated along 
PA 50°. The two bright lobes in the image are separated 
by 210 mas (30 AU) with the south-west side being signif- 
icantly brighter. A two-dimensional Gaussian fit to the 
visibility data gives a FWHM 244±6 mas x 123±9 mas 

^ All calibrators were assumed unresolved ex cept a Tau where w e 
adopted a uniform disk diameter of 19.8 mas (jPerrin et al.lllQQSh . 



along PA 49±1°. For a simplistic "flat disk" model, this 
2-to-l ratio elongation suggests we are viewing this disk 
oriented at ^60° inclination. 

A more suitable model for a circumbinary or transi- 
tional disk would be a tilted and asymmetric ring of emis- 
sion that is nearly unresolved along the minor axis (for 
example, see dramatic case of HR 4796; iKoerner et al.l 
119981: iSchneider et al.lll999D. Following the parameteri- 
zation of lMonnier et al.l (|2006l ) , we have fitted a "skewed 
asymmetric ring model" to our interferometry data here 
to better estimate the inner hole diameter. The best fit- 
ting model had an inner- hole diameter of 247x121 mas 
(35x17 AU) along PA 53°, with a 40% skew along PA 
284° (the thickness of the ring followed a gaussian profile 
with a FWHM of 25% ring radius). For purposes of 21 
we will assume the central hole in the circumbinary disk 
has a radius of 17.5 AU. 

Importantly, we want to emphasize that the scale of the 
emission is more than 4x larger than the separation of 
the binary at the heart of this system, and the elongation 
is along a distinctly-different position angle. Since our 
mid-IR data in 2004 were taken within one week of the 
near-IR data, we can clearly prove that the elongated 
mid-IR emission comes from a distinct and much larger 
component than the near-IR. We discuss the relationship 
between the circumbinary disk and the underlying binary 
orbit in gl 

At f irst glance, it may seem surprising that iLiu et aP 
(|2005l ) did not more clearly resolve this large circumbi- 
nary disk using their nulling technique, finding a FWHM 
of only 20 AU. However, this nulling measurement was 
performed only for PA 164°, which is along the narrow di- 
mension of the elongated emission we detected. Another 
possible explanation for the smaller size is that the band- 
pass filter used here has more contribution from 11.2/im 
PAH emission, which could be more extended than the 
mid-IR continuum. 

3.1.2. Near-infrared: Binary star 

Figure [2] shows the BSMEM image of v892 Tau at 
2.2/im. We confirm the binary nature of this target as 
first reported by S mith et al.l (|2005l ). We also set a limit 
of <10% of the emission possibly com ing from any sort of 
extended or "halo" component (e.g;., iKataza fc Maiharal 
119911 : iHaas et al.l 119971 : iLeinert et al.ll200lf ). In order to 
extract the separation and position angle of this binary, 
we have used both image fitting and direct fitting to the 
inter ferometric observables to yield the following result: 
p =44.2±1.0 mas, =79.9±1.0°, Flux Ratio 1.15±0.04 
(west component brighter than east component). 

3.2. Binary orbit 

We combined our new binary measurement (at 
2004. 67) with the 1996 .75 and 2003.76 measurements 
from 'S mith et al.l (|2005l ) and these data are plotted in 
Figure [3l The sm all proper motion (A ~10 mas= 
1.4 AU) observed by I Smith et al.l ([2005) over 7 years is 
problematic for this object since this would require a long 
period orbit (>100 yrs) inconsistent with the spectral 
types and luminosities of the binary stars themselves. 

Alternatively, the binary period could be approxi- 
mately 7 years, meaning the stars had gone nearly ex- 
actly once around their orbit between epochs. Although 
the timing is suspicious, the scenario is plausible except 
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this would require a system mass >20 Mq, much too high 
for the B8 spectral type and system luminosity. 

The scenario we favor is that infrared variability af- 
fected the brightness ratio between 1996 and 2003, con- 
fusing the assignment of "primary" and "secondary" star. 
This variation could be intrinsic or be due to varying 
line-of-si^ht obscu ration through this asymmetric disk. 
ISmith et al.l (|2005f ) did report a change in flux ratio, flnd- 
ing a brightness ratio of 1.0 for 2003.76 which implies that 
their reported position angle had a ±180° ambiguity. 

In our 2004.67 data, we definitively detect the south- 
west component as 15% brighter, consistent with the PA 
assignment of the published 2003.76 measurement. By 
flipping the position angle of the earlier 1996 measure- 
ments, we flnd a robust family of orbits with periods of 
~14 years, compatible with the measured binary sepa- 
rations and the expected system mass. Note that future 
orbital reflnement will lead to precise stellar masses for 
these intermediat e-mass, pre-main-seq uence stars (see 
also MWC 361 A. iMonnier et al.l[2006h - a rare oppor- 
tunity to advance the calibration of pre-main-sequence 
tracks for this mass range. 

To estimate the orbital parameters for v892 Tau, 
we used Monte Carlo sampling of the measured stel- 
lar separations along with a system mass constraint of 
5.5±0.5M(7) expected f or two stars of spectral type B8V 
(|Palla fc Stahled Il993f ). Figure [3] shows the main re- 
sults of our orbital study: P = 13.8 ±1.5 yrs, a = 
72.4 ± 6.3 mas (10 AU), e = 0.12 ± 0.05, i = 60.6 ± 3.8°, 
00 = 233 ± 42°, (^ = 28 ± 5°, To = 55480 ± 900 MJD. 

4. DISCUSSION 

We interpret the mid-IR emission of v892 Tau as a cir- 
cumbinary disk based on several arguments. The mid-IR 
emitting region is much larger than the separation be- 
tween the binary components, ruling out models which 
have all the mid-IR emission coming from disks around 
the indi vidual stars. In theoretical m odels of tidal trun- 
cation (|Artvmowicz fc Lubowl Il994l ) , the circumbinary 
disk has an inner hole approximately 1.8-2.6 x larger than 
the semi-major axis of the binary system (for eccentricity 
between and 0.25). For comparison, our best estimates 
of the hole radius and orbital semi-major axis are 17.5 
AU and 10 AU respectively - giving a ratio of 1.75 which 
is close to the theoretical expectation. Additionally, the 
position angle and derived inclination angle of the mid- 
IR emission (PA~50°, i ~60°) is similar to those derived 
from the orbit (PA~28°, i ~61°) and the mild asymme- 
try in the mid-IR emission suggests dynamical interac- 
tions between an eccentric binary and the surrounding 
disk (via resonances or disk warping). 

We further test the circumbinary disk hypothesis by 
analyzing the spectral energy distribution (SED). The 
SED shows a large m id-IR bump similar to that seen in 
"transitional" disks ("Calvet e t al.l l2002l : lEspai llat et al.' 
|2007). By reddening template Kurucz spectra for B8V 
stars, we can flt the near-IR emission with mostly pho- 
tospheric light (some near-IR emission from hot dust is 
allowed but not well constrained) assuming Ay = 10.95 
yielding a combined stellar luminosity of 400 L©, reason- 
able for two B8V stars. Our new proposed SED decom- 
position for v892 Tau can be found in Figure [31 showing 
that the rising mid-IR flux can be roughly characterized 
as a T~ 450 K blackbody, consistent with emission from 



the warm inner wall of the circumbinary disk. 

Using equation (1) in llsella et al.l (|2006l ). the predicted 
radius for the warm inner wall for the above stellar lu- 
minosity (assuming 0.25 jam spherical silicate grains) is 
18 AU - consistent with our derived inner hole radius 
of 17.5 AU based on imaging. Furthermore, we flnd 
that the scale height of this wall at 18 AU is ~ 1.8 AU 
(jDullemond et al1l200ll ) which can produce the observed 
line-of-sight Ay ^ 11 for the disk inclination of ~65^. 
This latter phenomena may not be widely appreciated - 
circumbinary (or transitional) disks have warm puffed- 
up inner walls (the rim scale height is a stronger func- 
tion of radius than temperature) at large radii that ef- 
fectively increases the possibility that central stars will 
be obscured or at least reddened as seen by distant ob- 
servers. The lack of scattering nebulosity in archival 
Hubble Space Telescope images independently suggests 
that the high line-of-sight Ay is likely from absorption 
by dust in the outer disk and not infalling envelope ma- 
terial. While a detailed model is beyond the scope of 
this Letter, a careful study will allow precise constraints 
of stellar luminosities, dust properties, and circumbinary 
structure and should be straight-forward with today's 3- 
D Mont e Carlo radiative transfer codes (e.g., TORUS, 
lHarriesll200a) . 

5. CONCLUSIONS 

We have discovered an extensive circumbinary disk 
around v892 Tau in the mid-infrared. We also inde- 
pendently conflrm the binary nature of the underlying 
stellar system and our new measurement allows us to flt 
an astrometric orbit, flnding a period ~14 years for sys- 
tem mass of ~6 M© . Our limited orbital phase coverage 
and some ambiguity in position angles allow only crude 
estimates and we strongly urge continued monitoring of 
this system using near-IR speckle, aperture masking, or 
adaptive optics. 

We have proposed a new SED decomposition with a 
line-of-sight extinction {Ay ~ 11) higher than previously 
thought, implying a system luminosity ~4OOL0. This 
result highlights the fact that circumbinary disks (and 
transitional disks) have much larger opening angles since 
the puffed-up inner wall causes enhanced extinction of 
the central stars even at inclinations of 60°. 

v892 Tau is another case where high-resolution imaging 
has motivated a fundamental shift in our understand- 
ing of an individual object. As spectral energy distri- 
butions are used to discover "transitional" disks impli- 
cating planet formation, we must be cognizant of the 
important role of binarity and the associated circumbi- 
nary disks that can mim ic signs of planet formation (e.g., 
[Ireland fc KrausI l2008l ) . Far from merely being spoil- 
ers to planet flnders, new circumbinary disks offer fresh 
laboratories for studying dynamical interactions between 
gas-rich disks and the massive orbiting bodies embedded 
within them as well as critical opportunities to calibrate 
pre- main sequence tracks. 
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Fig. 1. — This figure shows visibihty data (color scale) for v892 Tau and the specific Fourier coverage (dots). The left panel shows results 
for 10.7/im using the Keck segment tilting method while the right panel applies to the 2.2/im Keck aperture masking data. 
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Fig. 2. — This figure shows the image reconstructions of v892 Tau at 10.7/im (left panel) and 2.2/im (right panel) using the BSMEM 
software. At the bottom-left of each panel, we have included an estimate of the resolution of each image, corresponding here to 80 mas at 
10.7/im and 16 mas at 2.2/im. Contour levels are shown for the extended emission in the mid-IR spaced logarythmically in factors of 4: 
25%, 6.3%, and 1.6% of the peak. The scale bar applies for a distance of 140 pc. For orientation. East points left and North points up. 
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Fig. 3. — (left panel) This figure shows the separation and position angle measurements for v892 Tau from Smi th et all (|2Q05l ') and this 
work. We have plotted the best-fit model for a system mass of 5.5 M© (period 13.8 yrs, eccentricity 0.1 2) w ith the solid curve. Note 
that we had to fiip the position angle of the 1996 epoch measurement in order to find a realistic orbit (see ^3.21 for full d escripti on of our 
procedure and justification), (right panel) This figure shows the spectral energy distribution for v892 Tau (photometry from Strom & Strom 
|T994; Acke fc van den Ancker__20Q4 : Her bst & Shevchenko 1999, and our own JHK photomery from MDM Observatory) and our proposed 
decomposition of star plus disk. We have plotted the contribution from the reddened {Ay = 10.95, Ry = 6.4) stellar photospheres (2xB8V, 
system luminosity of 400 L©) and the contribution from the "warm inner wall" of the circumbinary disk with temperature T~450K. 



